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Ultra slow-roll inflation has recently been used to challenge the non-Gaussianity consistency re- 
lation. We show that this inflationary scenario belongs to a one parameter class of models and 
we study its properties and observational predictions. We demonstrate that the power spectrum 
remains scale-invariant and that the bi-spectrum is of the local type with / NL = 5(3 — n s )/4 which, 
indeed, represents a modification of the consistency relation. However, we also show that the sys- 
tem is unstable and suffers from many physical problems among which is the difficulty to correctly 
WMAP normalize the model. We conclude that ultra slow-roll inflation remains a very peculiar 
case, the physical relevance of which is probably not sufficient to call into question the validity of 
the consistency relation. 
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I. INTRODUCTION 

The theory of inflation convincingly describes the phys- 
ical conditions that prevailed in the very early Uni- 
verse [1-5]. However, there are many models of inflation 
and it is not yet clear which scenario is actually realized 
in Nature. For this reason, the recent developments in 
the calculations of higher correlation functions [6] are im- 
portant since they might allow us to constrain and maybe 
rule out many models of inflation. For instance, the sim- 
plest scenarios (i.e. a slowly rolling single field with a 
canonical kinetic term) are known to predict a negligi- 
ble level of non-Gaussianity, of the order of the slow-roll 
parameters [7-12]. Therefore, if any non-Gaussianity is 
detected in the future (for instance with the Planck satel- 
lite), these models would be excluded. 

Recently, however, it was argued in Ref. [13] that this 
is not necessarily true for the simple scenarios mentioned 
above. An explicit counter-example was investigated in 
Ref. [13] and it was shown that, in this particular case, 
the value of the / NL parameter can be a few instead of 
being negligible. Since this result challenges a well-known 
and important theorem, it is important to study in more 
detail the model that has been utilized to obtain this 
conclusion. In particular, one would like to know whether 
this just represents a very peculiar case or whether this 
can correspond to a generic class of meaningful models. 

In fact, the inflationary scenario used in Ref. [13] has 
been known for a long time and is named "ultra slow- 
roll" inflation. It was studied for the first time in Ref. [14] 
(Similar situations were also investigated in Ref. [15]). In 



the present article, we show that it belongs to a broader 
class of models that we explicitly identify. The goal of 
the paper is then to study this new family of scenarios, 
their properties and the corresponding observational pre- 
dictions (power spectrum and bi-spectrum). 

The paper is organized as follows. In the next sec- 
tion, Sec. II, we introduce ultra slow-roll inflation and 
show how it can be generalized. Then, we study the 
stability of the system and investigate whether one can 
easily produce 60 e-folds in the ultra slow-roll regime. 
Then, in Sec. Ill, we calculate the power spectrum of 
curvature fluctuations and show that it can be scale in- 
variant even if the slow-roll parameters are not all small. 
In Sec. IV, we estimate the non-Gaussianities and show 
that the Maldacena's consistency relation is indeed vio- 
lated. As a consequence the / NL parameter is of order 
one in this class of models. Finally, in Sec. V, we discuss 
in more details the difficulties of ultra slow-roll inflation 
and present our conclusions. 



II. ULTRA SLOW-ROLL INFLATION 

Let us consider an inflationary model with a single 
scalar field (with a standard kinetic term). The equa- 
tions of motion for (f> and for the Friedmann-Lemaitre- 
Robertson- Walker scale factor a(t) (t denotes the cosmic 
time and, in the following, a dot means a derivative with 
respect to t) are the Fricdmann and the Klein-Gordon 
equations, namely 



3Hcl + U 
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where H = a/ a is the Hubble parameter and Mpi the 
reduced Planck mass. The background evolution can also 
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FIG. 1: Left panel: new family of ultra slow-roll potentials for different values of the parameter n. Right panel: classical ultra 
slow-roll and slow-roll trajectories for n = —3.01 (solid green line and dotted blue line) and n = —2.99 (dashed pink line and 
dotted dashed red line). The initial condition for the scalar field is chosen to be 0i n i = 0.1 (f>n m . 



be characterized in terms of the slow-roll parameters (or 
horizon-flow parameters) defined by 



din e. 



(3) 



where TV denotes the number of e-folds, N = ln(a/ai n i) 
(aim being the scale factor at the beginning of inflation). 
The hierarchy starts with eo oc 1/H which implies that 
the first slow-roll parameter can be expressed as 



ei 



H 



H 2 2M 2 l H 2 ' 



(4) 



The second slow-roll parameter can be used to express 
the acceleration of the field, namely 



e 2 = 2ei +2-!-r. 

H4> 



(5) 



Inflation requires t\ < 1 and the slow-roll approximation 
is valid if all the horizon-flow parameters are small, e,; <C 
1 during inflation. 

As discussed in the introduction, if the potential is ex- 
actly fiat, then the Klein-Gordon equation implies that 
(j)/(H<t>) = —3 and this corresponds to the situation dis- 
cussed in Ref. [14] and named "ultra slow-roll inflation". 
In this case, despite the flatness of the potential, the slow- 
roll parameters are not all small: usually ei -C 1 but ob- 
viously £2 = O(l). As a consequence, one could expect 
the power spectrum to deviate from scale- invariance but, 
as shown in Rcf. [14], and as discussed in more detail 
in the next section, this is in fact not the case. This 
makes this model a priori interesting since this shows 



that scale invariance can be obtained even if the slow-roll 
approximation is violated. This also raises the question 
of whether this is peculiar to the property = or 
whether this can also be obtained in a broader context. 
In order to investigate this issue let us consider the more 
general condition 



nHcj), 



(6) 



where n is now an arbitrary number, not necessarily equal 
to —3. Let us also assume that the first slow-roll param- 
eter is still very small, t\ <C 1. Obviously, the case n ~ 
corresponds to slow-roll and n = —3 to ultra slow-roll. 
The corresponding equations of motion are given by 



~ V (</>), 
{n + ^H^ + V^ = 0. 



(7) 

(8) 



From these equations, it is easy to check that (j> oc a", 
which implies that 



ei oc a 



2n 



e 2 ~ 2n, e 3 = 0. 



(9) 



In particular, for n = —3, one recovers the well known 
scaling ei oc a~ 6 , see Ref. [14]. 

From the two equations of motion (7) and (8), it is also 
straightforward to integrate the classical trajectory. One 
obtains 



n + 3 r+ V_ 



(10) 



where 4>- m i denotes the initial value of the inflaton. How- 
ever, it is not obvious that this solution will satisfy the 
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FIG. 2: Left panel: exact (numerical) evolution of the field (solid green line) compared to the ultra slow-roll solution (dashed 
blue line) and to the slow-roll solution (dotted red line). The parameter n is taken to be n = —2.99 and the initial condition is 
the same as in Fig. I (right panel), namely cf>- m i — 0.1 <f>u m . The inset shows the global evolution of the system on a larger time 
scale. Right panel: same as left panel but for n — —3.01. 



condition (6). Requiring that this is the case, we find 
that the potential must obey the following differential 
equation 



V 2(ti + 3) V V 



Vcf, \ ' n(n + 3) 



0. (11) 



Interestingly enough, this differential equation can be in- 
tegrated and leads to the following potential 



V((f>) = AI 4 



'n(2n + 3) 



6 Mpi 



2(n + 3) 
2n + 3 



(12) 



where M is an arbitrary mass scale to be fixed by the 
Wilkinson Microwave Anisotropy Probe (WMAP) nor- 
malization and 0o an arbitrary constant that, without 
loss of generality, we can take to be (j>o = 0. The po- 
tentials in Eq. (12) represent a new family of model de- 
pending on one parameter, n. These potentials are rep- 
resented in Fig. 1. If n < —3, then they are defined 
only in the range —<f>u m < 4> < 4>n m with <fi\i m /Mp\ = 
7r/2y / 6/[n(3 + 2n)]. It is clear that if n ~ —3, the po- 
tential is extremely flat, justifying the name "ultra slow- 
roll". 

Using Eq. (10), one can compute the classical trajec- 
tory exactly. Inserting Eq. (12) into Eq. (10) leads to the 



following result 



<t> VSR (N) =Mp M 



G 



n(2n + 3) 



x arcsm 



e-sin J n(2n + 3) * 



6 M 



pi 



(13) 



One can check that <j> = implies TV = 0. Let us notice 
that this expression is not well-defined for the slow-roll 
case n = 0. In this situation, one should use the following 
expression 



</> SR (A0=M PM 



6 



n(2n + 3) 



x arcsm 



e n(„+3W3 sin 



'n(2n + 3) 4> h 



6 



Mpi 



(14) 



The ultra slow-roll and slow-roll trajectories are repre- 
sented in Fig. 1. The interpretation of these results can 
be easily understood. The slow-roll solutions just fol- 
low the curvature of the potential. Therefore, if n < —3 
then the vacuum expectation value of the field increases 
(the field escapes at infinity and will meet the singularity 
at (f) = 0ii m ) and inflation proceeds from the left to the 
right while, if n > —3, the field value decreases toward 
the minimum of the potential and inflation proceeds from 
the right to the left. The ultra slow-roll solutions behave 
in a different manner. Firstly, they are very similar what- 
ever the sign of n + 3 provided n ~ —3 and, secondly, the 
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FIG. 3: Top left panel: numerical (exact) evolution of the first horizon-flow parameter ei (solid green line) compared to its 
ultra-slow-roll behavior (dashed blue line) for n — —2.99 and <j!>ini = 0.1 0ii m . Top right panel: same as top left panel but for 
the choice n = —3.01. Bottom left panel: Evolution of the quantity 8 = (j>/(H(j>) for n = —2.99 and (f>- m i = 0.1 4>nm- Bottom 
right panel: same as bottom left panel but with n — —3.01. 



field always asymptotically approaches the minimum of 
the potential (i.e. <j> = 0). In the case n < —3, this means 
that the field actually climbs up the potential. This is 
of course due to the fact that, initially, it possesses a 
non-vanishing and non-negligible velocity. 

In order to investigate the stability of the system, we 
have also numerically integrated the exact equations of 
motion. Recall that the ultra slow-roll solution has been 
obtained from the exact equations of motion by neglect- 
ing the kinetic term in the right hand side of the Fried- 
mann equation. This term, although very small, repre- 



sents a perturbation for the ultra slow-roll solution. It 
is therefore interesting to study whether the system can 
stay in ultra slow-roll during a large number of c-folds. 
The result is presented in Fig. 2. After a few e-folds the 
exact solution (solid green line) leaves the ultra slow-roll 
solution (dashed blue line). On a larger time scale (see 
the insets in Fig. 2), we see that for n > —3 (left panel), 
the field passes through the minimum, becomes negative 
and starts to climb up the potential in the region <f> < 0. 
Then it reaches a maximum, turns back and decreases 
toward the minimum. Obviously, this evolution is very 
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FIG. 4: Number of e-folds at which the ultra slow-roll solution 
is left as a function of the initial value of the field. The exact 
numerical result (solid black line) is in excellent agreement 
with the analytical estimate of Eq. (21) (dashed red curve). 



different from the ultra slow-roll one. For n < —3 (see 
the right panel), the field also becomes negative but, since 
the curvature of the potential is now negative, it simply 
escapes to infinity in the region <j) < 0. 

It is also interesting to study the behavior of the first 
slow- roll parameter and of S = <fr/(H(j)). They are rep- 
resented in Fig. 3. The conclusions obtained before are 
confirmed. We see that e\ (top panels) scales as a 2n 
only for a few e-folds and then leaves the ultra slow-roll 
regime. For n = —2.99 (top left panel), we also notice 
that €\ vanishes and, of course, this corresponds to the 
point where, in the region (f> < 0, the field reaches a max- 
imum and turns back. This is confirmed by the fact that, 
in the case n = —3.01 (top right panel), the above men- 
tioned behavior never happens. Then, after this tran- 
sitory regime, in both cases, ei becomes constant with 
a very small value. The behavior of S (bottom panels) 
can be understood in a similar fashion. Initially S ~ —3 
since the field starts from the ultra slow-roll regime. Af- 
ter a few e-folds, this solution is left and, eventually, S 
reaches a regime where it remains constant with a very 
small value. In the case n = —2.99 (bottom left panel), S 
diverges when (f> — while its evolution remains smooth 
if n = —3.01. It is clear that having e± and 5 (or, equiva- 
lently £2) small and constant corresponds to nothing but 
the slow-roll regime. The conclusion of our numerical in- 
vestigation is therefore that the ultra slow-roll regime is 
unstable and is left after a few e-folds. Then, the system 
simply converges toward the slow-roll solution. 

It is also interesting to understand when the ultra slow- 
roll regime is left and what are the quantities which con- 



trol the instability. We now analyse this question in more 
detail. For this purpose let us define the following quan- 
tity 



nH(f> 



s_ 

5 

11 



(15) 



which is one during ultra slow-roll inflation. Using the 
equations of motion, it is easy to show that it obeys the 
following first order non linear differential equation 



dN 



<t><t> 



nH 2 



— + /(d - nf - 3) (16) 
n 



This equation cannot be solved exactly but we can study 
the behavior of small perturbations. For this reason, we 
now define A by mean of the following formula / = 1 + A. 
This quantity obeys the equation 

^ = £±2 Cl _ A [3 - e x + n(A + 2)] , (17) 
diV n 

which, in the regime where A -C 1, can be approximated 
by 



dA 
dN 



ei-A(3 + 2n) 



(18) 



Taking into account the behavior of the first slow-roll pa- 
rameter during the ultra slow-roll regime, namely ei = 
eiliniO 2 ", it is straightforward to obtain the following so- 
lution 

A(N) 



-(2ra+3)W 



n(4n + 3) 

For |n + 3 1 <C 1, one can approximate this solution by 

3 



(19) 



A(N)c- 



27 



(20) 




This allows us to estimate at which e-folds, Nd cv , the 
actual solution deviates from the ultra slow-roll one. 
Straightforward manipulations lead to 



(21) 



where (j>ini is the initial value of the field and A cr j an 
arbitrary value at which we estimate that one has left 
the ultra slow-roll solution. In the following, we estimate 
that this is the case if the actual solution differs for more 
than 10% from the ultra slow-roll one, that is to say 
A cr i ~ 0.1. We have computed this quantity numerically 
and have compared it with Eq. (21) in Fig. 4. Clearly the 
agreement is excellent. The main information brought 
by Eq. (21) is that the dependence in <^j n j is logarithmic. 
The ultra slow-roll solution is interesting if the system 
can follow the corresponding trajectory for at least 60 
c-folds. Using Eq. (21), one can estimate what it means 
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for the initial conditions. Straightforward manipulations 
lead to the constraint 



(22) 



In other words, in order to have 60 e-folds of ultra slow- 
roll inflation, one must fine-tune dramatically the initial 
value such that it is extremely close to the top of the 
potential. This is of course due to the logarithmic de- 
pendence in Eq. (21) which is in fact a consequence of 
the instability of the system. 



III. ULTRA SLOW-ROLL POWER SPECTRUM 

The fact that one of the slow-roll parameters is not 
small immediately raises the question as to whether the 
model can lead to an almost scale invariant power spec- 
trum. To address this question, it is convenient to 
work in terms of the so-called Mukhanov-Sasaki vari- 
able Vk, which is related to the curvature perturbation 
by Cfe = v kl {V2^Ip\(i\f£i)- The spectrum of can be 
expressed as 



2k 3 



8tt 2 M^ 



(23) 



The variable Vk obeys the equation of a parametric os- 
cillator, the time-dependent frequency being determined 
by the dynamics of the background [16] 



k 2 - — 

z 



i'k = 0, 



(24) 



where the prime denotes a derivative with respect to con- 
formal time and where z is given by a^JT{. The quantity k 
represents the comoving wavenumber of a Fourier mode. 
The "effective potential" z"/z can be expressed as 



z 
z 



a 2 H 2 



1 



1 



1 



ei 



: £ 2 



-£2^3 



(25) 



Despite the appearance of the slow-roll parameters, this 
expression is exact. As usual, the initial conditions on the 
perturbations are imposed when the modes are well inside 
the Hubble radius during inflation. In this regime, the 
modes do not feel spacetime curvature and, consequently, 
are usually chosen to be in the Bunch-Davies vacuum. 
This amounts to demanding that the Mukhanov-Sasaki 
variable Vk reduces to following Minkowski-like positive 
frequency mode in the sub-Hubblc limit: 



lim Vk = 

k/ (a H)— >oo 



1 



'2 k 



c 



— ikrj 



(26) 



In ultra slow- roll inflation, using Eq. (9), the effective 
potential for the perturbations can be expressed as 



1 



3n - 



(27) 
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FIG. 5: Spectral index versus parameter n for the new family 
of potentials. 



where r\ denotes the conformal time. We see that the 
solution to the mode equation can still be expressed as a 
Bessel function as it the case in the conventional situa- 
tion. The result reads 



Vk 



fa) = -M'e^i/^-M, (28) 



where H„ (z) is the Hankel function of first type. Then, 
for n < —3/2, the power spectrum on large Hubble scales 
can be written as 



V c (k) 



H 2 



7reiMp[ \aH 



k 



2n+6 



-77 ^usrW. 



where 



•^USFlfa) 



-2n-7 



T 2 (n + 5/2) cos 2 (n7r)' 



(29) 



(30) 



To our knowledge, this solution is new although the case 
n = —3 was found before in Rcf. [14]. If n > —3/2, then 
one has 



VAk) 



H 2 



7reiA/pj \aH 



-2n 



where 



9-l+2r 



T 2 (-n- 1/2) cos 2 (mr)' 
Finally, it remains the case n = —3/2. One finds 

^)^(^) 3 4> 2 (^ ^ 



(31) 



(32) 
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In all these expressions (and this is of course crucial for 
the case n < —3/2), ei must be evaluated not at the 
time of Hubble radius crossing but at the time of consid- 
eration, typically the end of inflation (of course, in the 
slow-roll case, this does not make a difference since the 
slow- roll parameters remain small and constant). The 
above expressions lead to the following spectral index for 
the power spectrum 

!2(n + 3), n < -3/2 

-2n, n > -3/2 

3 + 21n _1 [k/(aH)} , n = -3/2. 

The spectral index versus the parameter n is represented 
in Fig. 5. One sees that scale invariance is achieved for 
two values, namely n ~ which corresponds to the usual 
slow-roll and n ~ —3 which corresponds to ultra slow- 
roll. If n < —3 the spectrum is red while if n > —3, 
it is blue. It is easy to check that 0.96 < n s < 1, see 
Rcf. [17], corresponds to —3.02 < n < —3. Therefore, 
we obtain a new family of solutions leading to an almost 
scale invariant power spectrum but, clearly, n cannot de- 
viate from —3 too strongly. One can also re-express the 
spectral index in terms of the slow-roll parameters. For 
the slow-roll regime one obtains n s = 1 — ti while for the 
ultra slow-roll regime one has 

n s = 2?i + 6 = 7 + e 2 . (34) 

This should be compared to the standard slow-roll for- 
mula, n s — 1 — 2ei — £2. Of course, in the slow- roll 
regime, we obtain exactly the same equation given that 
ei -C 1- In the ultra slow-roll regime, however, we observe 
a breakdown of this result. This was already noticed in 
Ref. [14] for the case n = — 3 and it was shown in that 
reference that this is due to a breakdown of the horizon 
crossing formalism. Indeed, for n = —3, the slow- roll 
formalism leads to n s = 7 instead of the correct result 
n s = 1. 

It is also interesting to discuss in more detail the be- 
havior of curvature perturbations on large scales. Dur- 
ing inflation, the super-Rubble condition k/(aH) <C 1 
amounts to neglecting the k 2 term with respect to the 
effective potential z" / z in the differential equation (24). 
In such a case, it is straightforward to show that the 
super-Hubblc solution to v k can be written as follows 

v h (rj) ~ A k z(n) + B k z(n) J -^-y (35) 

where A k and B k are fc-dependent constants that are 
determined by the Bunch-Davics initial condition (26) 
chosen in the sub-Hubble limit. In our case, it is easy to 
show that this reduces to 

Cfc <xA k + B k a^ 2n+3 \ (36) 

In the slow-roll regime, the first term represents the grow- 
ing mode while the second one corresponds to the decay- 
ing one. In the ultra slow-roll regime however, the second 
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FIG. 6: Exact (numerical) evolution of curvature perturba- 
tions (solid blue line) versus the number of e-folds. The dotted 
dashed pink line represents the scaling oc a~'™ +1 ' while dot- 
ted dashed green line correspond to the scaling oc a~' 2n+3 \ 
The dotted red line is the quantity k 2 /(aH) 2 for a mode such 
that fc/ai n i ~ 50-Hi n i at the beginning of inflation. The dashed 
green line represents the quantity rfz"/z. 

term dominates over the first one (A similar situation was 
also studied in Refs. [18, 19]). This implies in particular 
that the power spectrum is still a time-dependent quan- 
tity on supcr-Hubble scales contrary to the standard case 
where it is conserved. This is apparent in Eq. (29) where 
the ei term in the denominator is a time-dependent quan- 
tity. On the contrary, the same factor in Eq. (31) is con- 
stant in time and, as a consequence, the slow-roll power 
spectrum does not evolve on large scales. It is also worth 
mentioning that curvature perturbations grow on sub- 
Hubble scales as well. Indeed since Q k ~ V k / (flv/ei) and 
since \v k \ stays constant in this case, this immediately 
implies \( k \ oc a~( n+1 \ In the slow-roll case, curvature 
perturbations decreases oc a . 

In order to check these considerations, we have nu- 
merically integrated Eq. (24). The result is presented in 
Fig. 6. The modulus of curvature perturbations corre- 
sponds to the solid blue line. The effective potential for 
the perturbations i] 2 z"/z is the dashed green line while 
k 2 /(a 2 H 2 ) is the dotted red line. When the dotted red 
line is above the dashed green one, the mode is within the 
Hubble radius and when it is below, the mode is outside 
the Hubble radius. In Fig. 6, we see that the mode starts 
its evolution deep inside the Hubble radius and crosses 
it out around N ~ 6. We verify that, inside the Hub- 
ble radius, \^ k \ grows like a- (n+1) , this particular scaling 
being represented by the dotted-dashed pink line. When 
the mode crosses out the Hubble radius, it is apparent 
that the behavior of \( k \ is modified. The dotted dashed 
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green line represents the scaling a _ ^ 2 ™ +3 ' and one sees in 
the figure that it is indeed the scaling of | Cfe I - Therefore, 
our numerical integration confirms that, in ultra slow- 
roll inflation, curvature perturbations grow on small and 
large scales. Around N ~ 13, ultra slow-roll inflation 
comes to an end and, as a consequence, the growth of 
Cfc stops. Then, as clearly seen in the figure, stays 
constant as usual in the slow-roll regime on large scales. 

This continuous growth of curvature perturbations 
during ultra slow-roll inflation turns out to have im- 
portant physical implications. Since V^(k) is a time- 
dependent quantity even on large scales, this means that 
the amplitude of the power spectrum at the time when 
inflation ends must now be compared with the WMAP 
normalization (in the slow-roll case, it is sufficient to nor- 
malize the power spectrum when the modes of cosmolog- 
ical interest today leaves the Hubble radius during infla- 
tion). If n ~ —3, the power spectrum of the curvature 
perturbation at the time when inflation ends is given by 

p < = ^ir 6a ' ¥ -(0- (37) 

where AJV, ~ 50 — 60 is the number of e-fold between 
the Hubble radius crossing time of the relevant mode 
and the end of inflation. The quantity ei* is t\ evaluated 
at the Hubble radius crossing time. From the WMAP 
normalization 7\- = 2.4 x 10 -9 , see Ref. [17], we find 
that for A A* = 60, M must satisfy 

which is far below the Big Bang Nucleosynthesis (BBN) 
bound M > O(McV). The result is expected. The quan- 
tity |£fc| grows so much during ultra slow- roll inflation 
that, in order to match the correct level of Cosmic Mi- 
crowave Background (CMB) fluctuations, one must com- 
pensate by a tiny mass scale in the potential. Let us 
notice that we also implicitly assume that, after infla- 
tion, the growth of stops. In addition, the above 



estimate is very conservative because it is expressed in 
terms of ei*. Since ei is decreasing from the beginning of 
inflation, it is likely that ei* -C 1. In other words, instead 
of AA*, the constraint could also be written in terms of 
the total number of e-folds. This means that a physically 
relevant ultra slow-roll inflation model can last only for a 
much shorter period than the 60 c-folds usually required. 



IV. ULTRA SLOW-ROLL NON GAUSSIANITY 

Let us now turn to the calculation of the three-point 
correlation function. For the case n = —3, the calculation 
was done for the first time in Ref. [13]. Here we general- 
ize this result for an arbitrary value of the parameter n. 
As is well-known, for slow-roll single field inflation with 
a standard kinetic term, the level of non-Gaussianity is 
very small, of the order of the slow-roll parameters, see 
Refs. [7-12]. This result is still true for ultra slow-roll 
inflation but, now, one of the slow-roll parameters is of 
order one. Therefore, one expects a / NL parameter of or- 
der one as well. We will see that this is what happened 
although, as noticed in Ref. [13], the relation between 
/ NL and n s is modified. 

The scalar bi-spectrum B s (fei, k 2 , ^3) is defined in 
terms of the three point correlation functions of the 
Fourier modes of the curvature perturbation £ as fol- 
lows [17, 20]: 

(CfeiCfc 2 Cfe 3 ) = (2rr) 3 6 s (fei,fe2,fc 3 ) 

x <5 (3) (fei + k 2 + ka) ■ (39) 

For convenience, we shall set G(ki,k 2 ,k 3 ) = 
(2 7r) 9 / 2 i? s (fci , k 2l k 3 ). Using the Maldaccna for- 
malism [6], the quantity G(ki,k 2 ,k 3 ) can be expressed 
as [21-23] (recall that the function below is the mode 
function that appears in front of the annihilation and 
creation operators in the canonical decomposition of the 
operator () 



G{k!,k 2 ,k 3 ) = ^2 G c (k 1} k 2 ,k 3 ) 



c=i 



C=l 



fkr (m) fk 2 (m) fk 3 (m)Q c ( fc i ,k 2 ,k 3 ) + (r? f ) (m ) /£ 3 (^ ) g* c (k x , fe 2 , fe 3 ) 



+GV(fci,fc 2 ,fc3), 



(40) 



where rjf denotes the final time when the bi-spectrum is to be evaluated. The quantities Q c (k\, k 2 , k 3 ) with C 
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1, • ■ • ,6 are described by the integrals [21-23] 



Gi(ki, fe 2 , k 3 ) 
G2{ki,k 2 ,k 3 ) 

G 3 (ki, k 2 , k 3 ) 
C? 4 (fei, fe 2 , fe 3 ) 
Gs{ki,k-2,k 3 ) 

Gefkx, fe 2 ,fc 3 ) 



rm 

2i di] a 2 e\ (/^ / fe * / fe * + two permutations) , 

-2i (ki ■ k 2 + two permutations) / dry a 2 e\ /£ a /£ 3 , 



rvt 

2i di] a 2 e\ 

J 71; 

rm 



kl k 2 2 ) fhi fk 2 f'k 3 + nvc permutations 



rm 

/ d?y a 2 ex e' 2 (/^ /£ 2 /£* + two permutations) , 
fei • fe 2 



d?7 cr ef 



/fei /fc* 2 /fe* 3 + nve permutations 



d?y a 2 el 



k\ (fe 2 ■ fe 3 



fc 2 k 2 



/fei fhi fks + tw0 permutations \ , 



(41) 
(42) 

(43) 
(44) 
(45) 
(46) 



where rji denotes the time when the modes are well in- 
side the Hubble radius during inflation. The additional, 
seventh term G7(fc 1 ,fc 2 ,fc 3 ) arises due to a field redefi- 
nition, and its contribution to G(fci, fe 2j ^3) is found to 
be 



immediately read / NL which is given by 1 



G r (fci,fc a ,fc 3 ) = f -2(2n + 3) 



two permutations 



\hAm)\ 2 \h 3 (m)\ 2 

(47) 



In the ultra slow-roll case, since e\ is very tiny while 
e 2 = 0(1), the above equations show that G7 gives the 
dominant contribution to the bi-spectrum for any config- 
uration of the triangle formed by the vectors k\, fc 2 and 
k 3 . Notice that the second term — 2(2n + 3) in Eq. (47) 
is absent in the standard slow-roll case. This originates 
from the fact that the terms in the cubic action that 
must be removed by field redefinition are of the form 
ae 2 C 2 /2 + 2((' /H + • • • , where the dots denote terms 
that always involve a spatial derivative of the curvature 
perturbation. In the standard case, only the first term is 
important because of the conservation of curvature per- 
turbations on supcr-Hubblc scales. On the other hand, 
in the present case where the decaying mode dominates 
over the growing mode, the second term also contributes 
since £' 7^ [13]. It is actually this second term that 
leads to the violation of the standard non-Gaussianity 
consistency relation. Then, the bi-spectrum becomes 



B s (k 1 ,k 2 ,k 3 



3 , (27Q- 1 / 2 
- (n + 2) ,3,3,3 



x [fcfP c (jfci)P c (*2) + 2 perms] (48) 



Interestingly enough, the bi-spectrum is of the local type 
not only in the squeezed limit but also for any other set 
of (fei, fc 2 , k 3 ). Then, from the above expression, we can 



f 1 ' 



SR 
NL 



(49) 



As noticed in Rcf. [13], this gives a relation between / NL 
and n different from the Maldacena consistency relation 
which yields = 5(1 - n s )/12 ~ 5n/6. Finally, it is 
also interesting to provide a relation between / NL and n s : 



/■ 



USR 
NL 



(50) 



where we emphasized again the fact that it is valid for 
any configuration, not only in the squeezed limit. This 
clearly shows that / NL becomes of order one even if the 
power-spectrum is almost scale invariant. Such a signal 
would be marginally detectable by the Planck satellite 
which, in principle, can see |/ NL | > 5. Finally, let also 
mention that, in order for the bi-spectrum we have just 
calculated to describe the non-Gaussianity which would 
actually be observed in the sky, it is necessary to assume 
that the growth of Cfc stops after the end of inflation 
and that reheating will not modify the result. The latter 
seems very reasonable as recently shown in Ref. [26] . 



V. DISCUSSION AND CONCLUSIONS 

Let us now recap our main results. Ultra slow-roll is 
not new and was studied in Ref. [14]. It is characterized 
by a situation where the first horizon flow parameter is 
very small but the second one is of order one. In this 
paper, we have generalized the ultra slow-roll regime to 



1 We use the same / NL as the one used by WMAP. Notice that 
Refs. [6, 22, 24, 25] use a different sign convention. 
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a one parameter family models. We have seen that, in ul- 
tra slow-roll inflation, the curvature perturbation can be 
dominated by the decaying mode. Despite this property, 
the corresponding power spectrum remains scale invari- 
ant and, hence, in agreement with the CMB observa- 
tions. This leads to the interesting situation where / NL 
is of order one even in a single field model with a stan- 
dard kinetic term. This clearly violates the Maldacena 
consistency relation. 

However, ultra slow-roll inflation appears to be 
plagued with many difficulties. Firstly, the system is un- 
stable and the ultra slow roll solution is left after a few 
c-folds only unless one artificially fine tunes the initial 
conditions. Secondly, the continuous growth of curvature 
perturbations implies that the mass scale of the potential 
must be extremely small in order to match the observed 
level of CMB anisotropy. In fact the corresponding value 
of M turns out to be unphysical. There is also a third dif- 
ficulty that we now discuss. As is well-known, when the 
potential is very flat, the quantum effects can dominate 
over the classical dynamics. In ultra slow-roll inflation, 
the typical variation of the scalar field (during one c-fold) 
due to the classical dynamics can be expressed as 



For A c 



0.1 and 



0.01, this gives 



A0 C 



-3A/& Vj 
3 + n V' 



(51) 



On the other hand, typical quantum jumps are given by 
A</>quant — H / (2tt) . Therefore, the classical equations 
of motion are valid only if A0 c j 3> A</> quant . Using the 
previous considerations, this leads to 



Mi 



> 



M 



pi 



27r|n|V3M| 1 ' 



(52) 



Given the requirement (22), one can have 60 e- folds of 
ultra slow-roll inflation in the classical regime only if 



Af 2 



27r|n|v / 3M| ] My |n + 3| 
that is to say 



M /648tt 2 A c 



M 



pi 



V |n + 3| 



1/4 



-45 



(53) 



(54) 



M < 1.1 GeV. 



(55) 



This is larger than the BBN bound M > O(MeV) but re- 
mains rather small. As we have seen the WMAP normal- 
ization provides a much tighter constraint on M. Never- 
theless, it is likely that in a realistic realization of ultra 
slow-roll inflation the quantum effects play a dominant 
role. 

It seems therefore difficult to produce 60 e-folds of in- 
flation in the ultra slow-roll regime. One can wonder 
whether the very flat region of the potential could only 
represents a limited part of the full potential. It seems 
however difficult to understand how the field could enter 
this part of the potential with the correct initial condi- 
tions 4> = nH(f>. Of course if = exactly, then the 
previous condition is true but this does not represents a 
realistic case as there will always be corrections, even if 
extremely small. In this case, moreover, the dynamics 
would be completely controlled by quantum effects. 

In conclusion, ultra slow-roll inflation represents an in- 
teresting playground but it remains a challenge to build 
a physically relevant model that would exhibit in this 
regime. In fact, this shows how robust the Maldacena 
consistency condition is. In order to violate it, we are 
forced to consider situations that appear to be plagued 
with many physical difficulties. 
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